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forces.  When  compared  to  the  form  of  the  exact  solutions, 
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Two  specific  forms  of  attrition-rate  coefficients  are 
considered,  allowing  for  different  maximum  effective  ranges 
of  the  two  opposing  weapon  systems.  These  coefficients  might 
be  used  to  model  a  constant-speed  attack  against  a  static 
defensive  position.  It  is  shown  that  for  these 
attrition-rate  coefficients,  the  Liouville-Green 
approximation  is  not  consistently  reliable  for  predicting 
force  levels,  and  yields  exact  results  only  under  certain 
restrictive  conditions.  ‘  Furthermore  it  was  found  that 
methodology  is  not  presently  available  to  accurately  predict 
from  Liouville's  normal  form  the  error  which  will  be 
incurred  by  invoking  the  approximation  in  a  specific 
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I. 


INTBODDCTION 


Although  it  is  well  recognized  that  combat  is  a  complex 
random  process,  a  model  that  provides  "closed-form” 
solutions  can  often  be  an  aid  to  the  military  analyst  by 
explictly  portraying  interrelationships  between  combat 
parameters.  Deterministic  differential  equations  are 
frequently  used  to  model  the  combat  process,  and  vary 
greatly  in  their  detail  and  level  of  complexity.  The 
military  analyst  must  be  able  to  recognize  significant 
variables  and  their  relationships.  An  idealized  model  is 
useful  for  giving  a  "first  cut"  portrayal  of  such 
relationships.  In  particular,  one  can  often  obtain  an 
analytical  solution  to  such  models.  An  analytical  solution 
that  can  be  expressed  in  terms  of  "elementary"  functions  can 
often  provide  insights  into  the  combat  process  and  identify 
factors  affecting  battle  outcome. 

The  use  of  differential  equations  to  model  warfare  was 
pioneered  by  F.  H.  Lanchester  in  1914  [Ref  8].  This  original 
work  was  quite  idealized  and  assumed  that  the  fire 
effectivenesses  were  constant  throughout  the  battle.  Many 
extensions  to  Lanchester 's  classic  formulation  have 
subsequently  been  developed,  to  include  such  factors  as 
mobility  of  forces  and  fire  effectiveness  that  varies  during 
the  battle.  Taylor  and  Brown  [Ref  13]  have  developed  a 
mathematical  theory  for  solving  Lanchester-type  equations 
with  these  temporal  variations  in  fire  effectiveness.  The 
fire  effectiveness  of  forces  is  represented  by  what  are  now 
called  Lanchester  attrition- rate  coefficients.  Taylor  and 
Brown  give  results  for  two  specific  forms  of  the 
coefficients:  (1)  effectiveness  of  each  side’s  fire 
proportional  to  a  power  of  time,  and  (2)  effectiveness  of 
each  side's  fire  linear  with  time  but  a  nonconstant  ratio  of 
attrition-  rate  coefficients. 
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Taylor  and  Brown  present  accurate  numerical  solutions 
to  equations  with  the  above  attrition- rate  coefficients. 
These  solutions  are,  in  general,  stated  in  terms  of  infinite 
series.  The  solution  method,  however,  is  rather  complicated 
and  involved  the  development  of  new  mathematical  functions 
alien  to  the  military  analyst.  Exact  results  are  certainly 
beneficial,  but  the  analyst  also  needs  a  solution  in  terms 
of  ”familiar''  functions  to  permit  parametric  analysis  and 
simplify  recognition  of  significant  variables.  To  provide 
such  a  solution  form,  Taylor  [Ref  14]  has  suggested  the  use 
of  the  so-called  Liouville-Green  approximation,  which  can  be 
expressed  in  terms  of  relatively  simple  functions.  For  an 
approximation  to  be  useful  in  predicting  force  levels, 
however,  it  must  be  reasonably  "close”  to  the  exact 
solution. 

This  paper  provides  a  numerical  evaluation  of  the 
Liouville-green  approximation  for  the  two  forms  of 
attrition-rate  coefficients  mentioned  above.  A  comparison  is 
made  between  the  approximation  in  each  of  these  cases  and 
the  numerical  results  given  by  Taylor  and  Brown.  Numerical 
examples  and  analytic  considerations  will  be  utilized  to 
investigate  the  accuracy  of  the  approximation  in  these 
cases. 


The  organization  of  this  thesis  is  as  follows.  We  first 
review  Lanchester's  classic  model  of  combat  between  two 
homogeneous  forces.  These  results  are  then  extended  to  the 
variable  attrition-rate  coefficient  case,  and  the  specific 
forms  to  be  studied  are  discussed.  The  representation  of 
the  solution  to  these  variable-coefficient  equations  in 
terms  of  so-called  general  Lanchester  functions  is 
presented.  Then  Liouville's  normal  form  and  approximation 
are  introduced,  evaluated,  and  discussed  for  the  two 
specific  forms  of  attrition-rate  coefficients. 
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II.  LANCHESTER ' S  CLASSIC  FORMOLAIION 


F.  S.  Lanchester  was  an  English  aeronautical  engineer 
who  first  developed  a  mathematical  formulation  for  combat 
between  two  conflicting  forces.  He  published  this  theory  in 

a  series  of  articles  in  the  British  journal  Engineering  in 

1914  [Ref  8j.  Lanchester 's  objectives  were  to  provide 
insights  into  "modern"  combat  and  to  quantitatively  justify 
the  military  principle  of  concentration.  He  stressed  that 
advances  in  technology  had  brought  long-range  delivery 
capability  of  modern  weapons  which  allowed  for  concentration 
of  firepower.  Lanchester  hypothesized  that  under  these 
modern  conditions,  combat  between  two  military  forces  could 
be  modelled  by: 

dx/dt  =  -ay  with  x  (t=0)  =  x  , 

0 

dy/dt  =  -bx  with  y  (t=0)  ~  r  (1) 

where  x(t)  and  y (t)  are  the  numbers  of  X  and  Y  forces  at 

time  t,  and  t=0  denotes  the  time  at  which  the  battle  begins. 

These  equations  are  valid  only  for  x  and  y  greater  than 
zero.  The  parameters  a  and  b  are  non-negative  constants  that 
have  come  to  be  referred  to  as  Lanchester  attrition-rate 
coefficients.  These  coefficients  represent  each  side's  fire 
effectiveness  or  firepower.  The  most  common  set  of 
circumstances  ,  under  which  these  equations  have  been 
hypothesized  to  apply  is  that  both  sides  use  aimed  fire  and 
target  acquisition  times  are  constant  [Ref  16]. 

From  (1)  Lanchester  deduced  his  classic  square  law 

2  2  2  2 
b(x^  -  X  (t))  =  a(y^  -  y  (t) )  ,  (2) 

which  has  the  important  implication  that  a  force  can 
significantly  reduce  its  casualties  by  initially  committing 
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more  forces  to  the  battle.  The  square  law  also  shows  that  in 
this  simple  model  the  Y  force  wins  (ie.,  X  is  annihilated) 
if  and  only  if 

X  /y  <  Va/b  .  ( 3 ) 
0  0 


From  (1)  Lanchester  also  derived  the  now  well-known 
results  for  the  time  history  of  the  force  levels 


X  (t)  =  X  cosh\/a1b  t  -  y  \Ja./h  sinhVab  t  , 
0  0 

y  (i)  -7  cosh  Vab  t  -  x  '^h/a.  sinhVab  t  , 
0  0 


(4) 


The  victory-predictio.n  condition  (3)  was  given  for  a 
"fight  to  the  finish".  As  H.  K.  Weiss  [Ref  15]  points  out, 
most  battles  terminate  before  complete  annihilation  of  one 
force  occurs.  Different  models  of  battle  termination  have 
been  proposed,  among  them  that  combat  ends  when  either  of 
two  given  "breakpoint"  force  ratios  is  reached.  Taylor  has 
shown  that  (3)  is  necessary  and  sufficient  for  I  to  "win"  a 
fixed-force-ratio-breakpoint  battle  [Ref  14]. 

Lanchester' s  original  formulation  (1)  involved  many 
implicit  assumptions  that  limit  its  usefulness  for  modelling 
real  combat.  The  forces  were  assumed  to  be  homogeneous, 
there  were  no  replacements  or  withdrawals,  and  no  movement 
of  forces  was  considered.  Setting  the  attrition-rate 
coefficients  a  and  b  in  (1)  as  constants  implies  that  the 
fire  effectiveness  per  firer  of  each  side  is  constant 
throughout  the  battle. 


III.  VARIABLE  ATTRITION-RATE  COEFFICIENTS 

When  modelling  mobile  weapon  systems  and  movement  of 
forces,  one  recognizes  that  the  assumption  of  constant  fire 
effectiveness  is  open  to  question  [fief  2].  If  fire 
effectiveness  is  permitted  to  vary  during  the  battle, 
Lanchester’s  formulation  (1)  takes  the  form: 

dx/dt  =  -a(t)  y  with  x  (t=0)  =  x  , 

0 

dy/dt  =  -b(t)  X  with  y  (t=0)  =  y^  ,  (5) 

where  a(t)  and  b  (t)  are  time-dependent  attrition-rate 
coefficients.  The  form  of  the  attrition-rate  coefficients 
depends  on  such  variables  as  force  separation  (range  of 
battle) ,  tactical  posture  of  the  targets,  firing  rate  and 
rate  of  target  acquisition  [Ref  5]. 

New  operations  research  techniques  for  predicting  these 
attrition-rate  coefficients  were  developed  in  the  1960 's.  A 
significant  contribution  in  this  area  was  the  development  of 
a  methodology  by  S.  Bonder  [Ref  2  and  4]  for  prediction  of 
coefficients  based  on  weapon  system  performance  data. 
Another  important  development  was  the  method  of  G.  Clark  for 
the  maximum  likelihood  estimation  of  these  coefficients 
based  on  Bonte  Carlo  simulation  output  data  [Ref  6].  The 
work  of  these  and  others  has  facilitated  the  use  of  models 
such  as  (5)  in  current  combat  studies. 

To  effectively  model  today's  highly  mobile  warfare,  the 
form  of  the  attrition- rate  coefficients  in  (5)  must 
accomodate  current  tactics  and  battle  conditions. 
Specifically,  the  coefficients  must  permit  the  fire 
effectiveness  of  one  or  both  sides  to  be  non-zero  at  the 
start  of  battle  (to  model,  for  example,  an  ambush).  The 
coefficients  must  also  allow  for  the  opposing  weapon  systems 
to  have  different  maximum  effective  ranges.  Since  mobile 
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forces  close  with  each  other  in  combat,  their  fire 
effectiveness  should  be  permitted  to  increase  or  decrease  as 
a  function  of  range  or  time.  The  attrition-rate  coefficients 
examined  in  this  paper  have  the  flexibility  to  model  these 
various  situations.  Motivated  by  previous  work  by  Bonder 
and  Farrell  [Ref  5],  Taylor  and  Brown  have  considered  the 
following  attrition-rate  coefficients: 

m  n 

a(t)  =  k  (t+C)  and  b(t)  =  k  (t+C+A)  .  (6) 

a  b 

The  specific  cases  of  (6)  to  be  considered  here  will  be  that 
in  which  A=0,  C>0, 

m  n 

a(t)  =  k  (t+C)  and  b(t)  =  k  (t+C)  ,  (j\ 

a  b  ^ 

and  that  in  which  m=n=1r  A>0,  C>0, 

a(t)  =  k  (t+C)  and  b  (t)  =  k  (t+C+A).  (8) 

a  b 

The  parameter  A  will  be  referred  to  as  the  ’’offset 

parameter”  since  it  permits  modelling  of  combat  between 

forces  whose  weapon  systems  have  different  maximum  effective 
ranges.  The  parameter  A  also  allows  the  ratio  of 
attrition-rate  coefficients  to  be  non-constant.  Koopman  [see 
Ref  10]  first  observed  that  if  the  ratio  of  attrition-rate 
coefficients  is  constant,  the  time  solution  to  (5)  is  no 
more  complicated  than  the  constant-coefficient  solution. 
That  is,  when 

a(t)  =  k  h(t)  and  b  (t)  =  k  h  (t)  ,  (9) 

a  b 

where  h (t)  is  the  common  time-dependent  factor  in  both 
attrition-rate  coefficients,  the  solution  is 

X  (t)  =  X  cosh  W  (t)  -  y  Jk  /k  sinh  W(t), 

0  0  a  b 
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where  W  (t)  =  \/  Ic  3c  •^h(s)  ds.  Koopman  also  found  that  the 

a  b  0 

constant  ratio  of  coefficients  again  yields  a  square-law 
relationship  (see  also  Ref  16). 

The  parameter  C  in  (7)  and  (8)  allows  for  the  initial 
fire  effectiveness  of  each  side  to  be  non-zero.  This  permits 
modelling  of  battles  which  begin  within  the  maximum 
effective  ranges  of  the  opposing  weapon  systems.  The 
parameters  m  and  n  allow  the  fire  effectiveness  of  each  side 
to  increase  or  decrease  with  time.  With  m,n<0,  the  fire 
effectivenesses  decrease.  This  situation  might  model  a 
battle  in  which  combatants  take  cover  and  improve  their 
positions  as  the  battle  progresses.  With  m,n>0,  the 
effectivenesses  increase,  such  as  might  be  the  case  in  which 
two  forces  are  closing  with  each  other. 

To  model  an  engagement  in  which  forces  close  with  one 
another,  it  is  beneficial  to  relate  their  instantaneous  fire 
effectiveness  to  the  range  separating  them.  As  this  range 
decreases,  one  might  expect  fire  effectiveness  to  increase. 
The  attrition-rate  coefficients  (6)  can  be  expressed  in 
terms  of  range  between  forces  and  the  maximum  effective 
ranges  of  the  opposing  weapon  systems: 

dx/dt  =  -a(r)  y  =  -a  (1-r/2  )  y, 

0  a 

dy/dt  =  "b  (r)  x  =  -b  (1-r/R  )  x.  (10) 

0  b 

These  equations  have  been  used  by  S.  Bonder  [Refs  1  and  3] 
to  model  a  constant-speed  attack  on  a  static  defensive 
position.  Figures  1  and  2  show  sample  relationships  between 
the  attrition-rate  coefficients  and  range  between  opposing 
forces  for  various  input  parameters.  In  these  equations,  r 
is  the  range  between  the  forces,  and  R  and  R  are  the 
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maximum  effective  raages  of  the  Y  and  X  weapon  systems. 


respectively.  The  constants  a 


0  0 


and  b  represent  the  fire 


effectiveness  of  each  side  at  r=0.  Range  and  time  are 

related  by  r(t)=R  -vt,  where  R  denotes  the  opening  range  of 

0  0 

battle,  and  v  is  the  constant  attack  speed.  With  this 
relationship,  the  attrition-rate  coefficients  in  (10)  can  be 
related  to  (7)  and  (8)  .  For  (7)  we  have; 


m  n 


k  =  a^(v/R  )  ,  k^  =  B  (v/R  )  ,  c 

a  0  a  b  0  b 


where  R  =  R  ^  R  ,  while  for  (8)  we  have: 
a  b  0 


k  =  a  v/R  ,  k  =  B  v/R  ,  A  =  (R  -R  )/v,  C  =  (R  -R  )/v, 
aOabOB  ba  aO 


where  R  >  R  >  R  . 
b  a  0 
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IV.  A  STANDARD  FORM  FOR  THE  VARIABLE-COEFFICIENT  SOLOTION 


The  variable  attrition-rate  coefficient  equations  (5) 
yield  the  force  level  equations 

2  2 

d  x/dt  -  C1/a(t)  da/dt  ]  dx/dt  -  a(t)b(t)  x  =  0,  (13) 
with  initial  conditions 


X  (t=0) 


X  and  r  Va  (t)  dx/dt]  =  -y  , 
0  t=0  0 


and  similarly 

2  2 

d  y/dt  -  [1/b(t)  db/dt]  dy/dt  -  a(t)b(t)  y  =  0,  (14) 

with  initial  conditions 

y  (t=0)  =  y  and  [1/b(t)  dy/dt]  =  -x  . 

0  t=0  0 

The  solution  to  (13)  is  given  by 

X  (t)  =  C  X  (t)  +  C  X  (t)  ,  (15) 

11  2  2 

where  [X  (t) ,  X  (t) ]  denotes  a  fundamental  system  of 
1  2 

solutions  to  (13)  ,  and  C  and  C  are  constants  determined  by 

1  2 

the  initial  conditions  [see,  for  example.  Ref  7].  The 

Y-force  level  has  similar  form.  The  functions  X  (t)  ,  X  (t)  , 

1  2 

Y  (t)  ,  and  Y  (t)  are  referred  to  by  Taylor  and  Brown  as 
1  2 

General  Lanchester  Functions.  Their  properties  are 

summarized  in  Table  I.  The  functions  X  (t)  and  Y  (t)  are 

1  1 

similar  to  the  hyperbolic  cosine,  while  X  (t)  and  Y  (t)  are 

2  2 

similar  to  the  hyperbolic  sine.  The  General  Lanchester 

Functions  may  be  constructed  by  successive  approximations  or 
by  infinite  series  methods. 
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A.  APPLICATION  TO  POWER  ATTRITION-RATE  COEFFICIENTS 


When  this  solution  method  is  applied  to  the  case  of 
power  attrition-rate  coefficients  (7)  ,  the  solution  (15)  maj 
be  written  as 


x(t)  =  X  [0  (C)  u  (t+C)  -  V  (C)v  (t+C)  ]  - 

0  m,n  m,n  m,n  m,n 

y  JTTiT  [u  (C)v  (t+C)  -u  V  (C)  ],  (16) 

0  a  b  m,n  m,n  m,n  m,n 


where  u  ,  v  ,0  and  V  are  referred  to  by  Taylor  as 
m«n  m,n  m,n  m,n 

Power  Lanchester  Functions.  They  have  the  properties  shown 
in  Table  II.  For  computational  convenience,  utilization  is 
made  of  the  Lanchester-Clif f ord-Schlaf li  (LCS)  functions 


00  2k 

F  (X)  =  Z  (x/2)  /  [k  n  (j+v)  ], 

V  k=0  j=0 

„  2k+1  k 

G  (X)  =  ?  (X/2)  /  [k  n  (j  +  v)  ]. 

^  k=0  jsO 

The  Power  Lanchester  Functions  may  be  expressed  in  terms  of 
the  LCS  functions  through  the  following  relationships: 


u  (t) 
m,n 


V  (t) 

m  ,n 

U  it) 
m,  n 


V  (t) 
m,n 


/  m  n 

where  S  (t)  W  (k  t  )  (k  t  )  t  / 


functions,  (16)  may  now  be  written  as 


F  (S(t)) 

q 

(m-n)/2 

t  G  (S(t)) 

P 

F  (S(t)) 

P 

(m-n)  /2 

t  G  (S(t)), 

q 


((m+n+2)/2).  Osing  the  LCS 
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X(t)  =  (J(0))P  (J(t))  -  (1+t/C)  (J(0))G  (J(t)) 

0  p  g  q  P 

-  y  Ja  (t=0)/b  (t=0)  [  (1+t/C)  (J(0))G  (  J  (t)  )  - 

0  q  p 


G  (J(0))F  (J(t)  )  ]. 

p  q 


(17) 


B.  APPLICATION  TO  OFFSET  LINEAR  COEFFICIENTS 


Similar  techniques  are  utilized  to  apply  this  general 
theory  to  the  case  of  the  offset  linear  attrition-rate 
coefficients  (8)  .  Hith  these  coefficients,  a  fundamental 
system  of  solutions  to  (13)  is  given  by 


x^(t)=f  (t+C)  , 


X  (t)=g{t+C),  y  (t)=F(t+C),  y  (t)=G(t+C) 
2  1  2 


where 


OO 

(^k 

2n 

(2n)! 

"  k  k  4n-k 

f 

(t) 

=  z 

{ 

k 

u/2) 

/ 

}  Z  B  A  t 

n- 

0 

a 

b 

k=0  n 

OO 

(^k 

2n+ 

1 

n  k  k  4n+2-k 

g 

(t) 

=  z 

C 

k 

/(2n 

+  1)  !  }  s  C  A  t 

n- 

0 

a 

b 

k=0 

OO 

(\/  k 

2n 

k  k  4n-k 

F 

(t) 

=  z 

{ 

k 

/ 

(2n)  ! 

}  ZD  A  t 

n= 

0 

a 

b 

k=0  ° 

OO 

2n+ 

1 

n+Jk  k  4n+2-k 

G 

(t) 

=  z 

[ 

(Vk 

k 

/(2n 

+  1)  !  }  Z  E  A  t 

n= 

0 

a 

b 

k=0  “ 

The  coefficients 


k 

B  , 

n 


k 

C  , 
n 


k 

D  , 

n 


k 

and  E  satisfy  complicated 
n 


recurrence  relations  [see  Taylor  and  Brown  (Ref  13) ],  The 

functions  f  (t)  ,  g  (t)  ,  F  (t)  ,  and  G  (t)  are  referred  to  as 
Offset  Linear  Lanchester  Functions.  As  in  the  case  of  the 
power  attrition-rate  coefficients,  use  is  made  of  the 
following  Auxiliary  Offset  Linear  Lanchester  Functions  for 
computational  convenience: 
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h  (p,u) 


oo  2n  ^  k  k 

=  Z  {p  /{2n)  !}  2  B  p 
n=0  k=0n 

00  2n+1  n  k  )c 

w(p,u)  =  S  CP  /(2n+1);}  ECU 

n=0  k=0^ 

00  2n  k  k 

H(p,u)  =  2  {p  /(2n)  !}  2  D  u 

n=0  k=on 

oo  2Q+t  n+ljc  k 

w  (PrU)  =  2  (p  /(2n+1) !  }  E  E  u  . 

n=0  k=0  ^ 

Otilizing  these  auxiliary  Lauchester  Functions,  the  X-force 
level  solution  (15)  may  be  written  as 

x(t)  =  X  [H  (J(0)  ,j  (0)  )  h(J(t),j(t))  -  W(J(0)  ,j  (0)  ) 

0 

‘'(J(t)  ,  j(t))  ] 

-  7jy^  /k  [h(J(0)  ,j  (0))  w(J(t),j(t))  -  W(J(0)  ,j(0)) 

0  a  b 

lx(J(t),j(t))  ],  (18) 

where  j  (t)  =  A/(t+C)  and  J(t)  =\Jk  k  (t+C)  /2. 

a  b 

The  equations  (17)  and  (18)  are  solutions  to  the 
X-force  level  equation  (13)  for  the  case  of  power  and  offset 
linear  attrition- rate  coefficients,  respectively.  This 
solution  method  of  Taylor  and  Brown  represents  a 
considerable  extension  of  Lanchester  theory  in  the  area  of 
variable  attrition-rate  coefficients.  Previous  analytic 
results  had  been  available  only  in  the  more  limited  case  of 
the  opening  range  of  battle  equal  to  the  minimum  of  the 
maximum  effective  ranges  of  the  two  weapon  systems.  That  is, 
the  initial  fire  effectiveness  of  at  least  one  side  was 
restricted  to  be  zero  [see  for  example  Ref  12]. 

While  the  time  solutions  (17)  and  (18)  appear 
formidable,  their  implementation  is  simplified  if  a  digital 
computer  is  available.  As  noted  earlier  in  this  paper. 
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however,  these  solutions  are  complex  and  involve  new 
mathematical  functions.  This  makes  it  difficult  for  the 
military  analyst  to  perceive  significant  relationships  among 
battle  variables,  and  renders  parametric  analysis  difficult 
and  time-consuming.  This  provides  motivation  for  the 
approximation  technique  to  follow. 
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7. 


REDUCTION  TO  L10aVILLE*S  NORMAL  FORM  AND  APPROXIMATION 


The  argument  ab  (referred  to  by,  among  others,  Taylor 
and  Brown  [Ref  13]  as  the  "intensity  of  combat")  in  the 
constant  coefficient  result  (4)  provides  motivation  for  the 
transformation  [Ref  14] 

t _ 

z  =  f\la.  (s)b  (s)  ds,  (19) 

^0 

where  z  (t=0)  is  denoted  z  .  Applying  this  transformation  to 

0 

(13)  and  (14)  yields 
2  2 

d  x/dz  +  CV2  d/dz  ln[b(t)/a(t)  ])  dx/dz  -  x  =  0,  (20) 
with  initial  conditions 


x(z=Zo) 


1/2 

X  and  C[b(t)/a(t)]  dx/dz}  = 

0  z=z  0 


and 


2  2 

d  y/dz  +  [1/2  d/dz  ln[a(t)/b(t)  ]}  dy/dz  -  y  =  0,  (21) 
with  initial  conditions 

1/2 

y(z=z  )  =  y  and  {[  a  (t) /b  (t)  ]  dy/dz)  =  -x  . 

0  0  z=z  0 


Taylor  and  Brown  have  shown  that  the  X  force  level 
equation  (13)  may  be  transformed  into  a  linear  second  order 
differential  equation  with  constant  coefficients  if  and  only 
if  (20)  is  a  constant  coefficient  equation.  This  is  the  case 
if  and  only  if 

( l/Va  (t)  b  (t) )  d/dt  ln[a(t)/b  (t)  ]  =  CONSTANT. 

Equation  (20)  shows  the  significance  of  the  parameters 
a(t)/b(t)  (the  relative  effectiveness  of  the  two  weapon 
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systems)  ,  and  the  intensity  of  combat  \/a  (t)  b  (t)  in 
determining  battle  outcome.  This  is  a  generalization  of  the 
well-known  constant  coefficient  results  to  the  case  of 
temporal  variations  in  fire  effectiveness.  The  significance 
of  these  parameters  may  be  seen  more  explictly  by 
transformation  to  Liouville's  normal  form  (see  Ref  14). 


The  transformation 

1/4 

x(z)  =  X(z)[  (a(t)/a^)  /  (b(t)/b^)]  ,  (22) 

may  be  applied  to  (20)  and  yields  the  normal  form  with  the 
first  derivative  of  the  dependent  variable  removed 


2  2 

d  X/dz  -  [  1  +  F  (z)  ]X  =  0, 


(23) 


with  initial  conditions 


X (z=z  )  =  X  and  dX/dz(z=z  )  =  -y  Va  /b  -  x  e  / 
00  0  00000 


where 


I  f 


P(z)  =  [R(t)  ] 


-1/4 


(24) 


F(z)  =  P  (z)  /  P(z) 

H(t)  =  a(t)  /b(t)  e(t)  =  1/4  \/a(t)  b  (t)  d/dt  InR  (t)  . 

(25) 

2  2 

In  (24),  the  notation  P*'(z)  denotes  d  P/dz  .  Equation  (23) 
is  Liouville's  normal  form  [see  page  23  of  Ref  9]. 


For  the  Y  force  level,  the  transformation 

1/4 

y(z)  =  Y(z)C  (b(t)/b^)  /  (a(t)/a^)] 
yields  the  normal  form 


2  2 

d  Y/dz  -  [1  +  G  (z)  ]Y  =  0, 


(26) 


with  initial  conditions 
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Y  (z=2  )  =  7  and  dY/d2(z=z  )  =  -x  /a.  +  y  e  , 

00  0  00000 


where 


1/4 

G(z)  =  Q"  (z)  /  Q(z)  and  Q  (z)  =  [  R  (t)  ]  =  1/P  (z) 


Expressing  the  transformed  force-level  equation  (23)  in 
the  simple  form 


2  2 

d  X/dz  -  X  =  Q  (z)  X, 

variation  of  parameters  may  be  used  to  obtain  the  solution 
to  (23)  as 

X  (z)  =  X  cosh(z-z  )  -  [y  \/a  /b  +  x  e  ]  sinh  (z-z  ) 

0  ^  0  00000  0 

+  /F  ( j)  sinh  (z- j)  X  ( j)  d  j .  (27) 

^0 

This  equation  may  be  expressed  in  terms  of  the  original  time 
variable  t  and  dependent  variable  x  as 

1/4  _ 

x(t)  =  [  (a  (t) /a  ) /(b(t) /b  ]  {x  cosh  (  (s)  b  (s)  ds) 

0  0  ^0  0 

-[  y  a  /b  +  X  e  ]sinh  (  /i/a  (s)  b  (s)  ds) } 

0  0  0  0  0  0 

1/4  1/4 

+  [a(t)/b(t)]  d/^s£1/(a(s)b(s)  )d/ds[b(s)/a(s)  3  } 

sinh  (  /i/a  (s  j)  b  ( j)  d  j)  X  (s)  ds.  (28) 

s 

If  the  relative  effectiveness  varies  slowly,  then  the 
relationship  (24)  would  indicate  that  F(z)<<1  so  that  one 
could  drop  the  integral  term  in  (27)  yielding 

X(z)  =  X  cosh  (z-z  )  -  [y  \I a  /b  +  x  e  ]  sinh  (z-z  ), 

0  0  00000  O^^gj 

which  is  the  Liouville-Green  approximation  [Ref  11].  The 
approximation  may  be  expressed  in  terms  of  the  original 
variables  t  and  x  as 
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X(t) 


1/4  ^  - - 

[  (a(t)/a  )/b(t)/b  ]  (x  cosh  (  /Va  (s)  b  (s)  ds) 

0  0  ^0  0 

-  [y  Va  /b  +  X  e  ]sinh  (  /\/a  (s)  b  (s)  ds) }  (30) 

0  0  0  0  0  0 
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VI 


APPLICATION  TO  POWER  OFFSET  LINEAR  COEFFICIENTS 


The  Liouville-Green  approximation  may  now  be  applied  to 
the  constant-speed  attack  model  with  attrition- rate 
coefficients  (7).  For  the  reader's  convenience,  these 
coefficients  are 


m 

a(t)  =  k  (t+C) 
a 


and 


b  (t)  =?  k  (t+C) 
b 


(7) 


This  model  assumes  that  the  maximum  effective  ranges  of  the 
two  opposing  weapon  systems  are  equal  (no  offset) .  The 
opening  range  of  battle  is  taken  to  be  within  this  maximum 
effective  range  (i.e.,  C>0) .  Applying  the  Liouville-Green 
approximation  (29)  ,  the  time  solution  to  the  X  force  level 
equation  in  terms  of  the  original  time  variable  t  and 
dependent  variable  X  becomes 


(m'-n)/4  - 

x(t)  =  (1+t/c)  (X  cosh/va  (s)  b  (s)  ds 

0  0 

, _  (m-n) /2  -(m+n+2)/2 

-  [y  Vk  /k  C  +(x  (m-n)/4  k  k  )C  ] 

0  a  b  ,.0  a  b 

sinh  /\/a  (s)  b  (s)  ds} , 

0 


(31) 


where 


/\/a  (s)  b (s)  ds  =  (Vk  k  /(m+n+2)/2) 
0  a  b 

(m+n+2)/2  (m+n+2)/2 

{  (t+C)  -C  }  . 


For  future  reference,  equation  (24)  becomes  in  this 

case 

2  2 

F  (z)  =  [  (m-n)  (3m+n+4)  ]  /  [4(m  +  n  +  2)  z  ],  (32) 

or  in  terms  of  the  original  time  variable  t  and  range 
parameters  (11): 
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m  n 

F(t)  =  [B  R  (m-n)  (3m+n+a)  ]  / 
a  b 

m+n  m+n+2 

[16a  b  (vei)  {t+(R  -H  )/vel)  ],  (33) 

0  0  a  0 


where  vel  =  constant  velocity  of  attack. 


He  can  also  apply  the  Liouville-Green  approximation  to 
the  case  of  the  offset  linear  attrition-rate  coefficients 


a(t)  =  k  (t+C)  and  b  (t)  =  k  (t+C+&)  .  (8) 

a  b 


In  this  case,  the  time  solution  to  the  X  force  level 
equation  beccmes 


1/4  - 

x(t)  =  [  (1+ (A/C)/(1+ (A/ (t+C) ) )  ]  (x  cosh/\/a(s)b(£)  ds 

0  0 

-  C  (jjk  /k  )/(!+  (A/C))  +  (X  (A/C))/ 

0  a  b  0 

2  3/2  t _ 

(4  k  k  C  (1+(A/C))  )  Jsinh/Va  (s)  b  (s)  ds}  ,  (34) 

a  b  0 


where 


/\/a  (s)  b  (s)  ds  =  (A  k  k  /8)  C  (1+ (2  (t+C)/A)  )  (1+ (2  (t +C) /A) )  -1 

0  a  b 

2 

-  ln[  1+ (2  (t+C) /A)  +  (1+ (2  (t+C)/A) )  -1]} 


For  future  reference,  equation  (24)  becomes  in  this 

case 


F(t)  =  R  B  (R  -R  )C12(t+(H  -R  ) /vel)  +  7  {  (R  -R  )/vel)]  / 
aBBa  aO  Ba 

3  3  3  ,  ^ 

{16a  b  vel  (t+ (R  -R  )/vel)  (t+ (R  -R  )/vel)  }.  (35) 

0  0  a  0  b  0 
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7II.  DISCaSSION 


A  computer  program  was  developed  in  FORTRAN  to  produce 
plots  of  the  X  force  level  versus  range  between  opposing 
forces  for  the  constant-speed  attack  model  described 
earlier.  Both  the  power  and  offset  linear  attrition-rate 
coefficient  cases  (equations  (31)  and  (34))  were  considered. 
These  approximate  solutions  were  compared  to  the  exact 
results  given  by  Taylor  and  Brown  [Ref  13]  to  evaluate  the 
accuracy  of  the  Liouville-Green  approximation  to  equations 
with  these  coefficients. 


A.  POWER  ATTRITION-RATE  COEFFICIENTS 

For  the  case  of  power  attrition-rate  coefficients,  the 

maximum  effective  ranges  of  the  two  weapon  systems  were  set 

equal  to  2000  meters,  with  an  opening  range  of  battle  of 

1250  meters.  The  fire  effectiveness  of  the  two  sides  at  zero 

range  separation  was  a  =.06  X  casualties/min-Y  unit,  and 

0 

b  =.6  Y  casual ties/min-X  unit.  The  velocity  of  attack  was 
0 

constant  at  five  miles  per  hour,  and  the  initial  force 
levels  were  X  =  10  and  Y  =  30.  These  parameter  values 

correspond  to  those  used  by  Taylor  and  Brown  in  their  paper 
[Ref  13].  The  exponents  m  and  n  in  (7)  were  varied 
throughout  the  range  explored  in  Ref  13.  The  results  are 
shown  in  Figures  3  through  10,  where  each  figure  contains 
the  approximate  time  solution  and  the  corresponding  exact 
analytical  solution  obtained  by  Taylor  and  Brown.  As  even  a 
cursory  examination  reveals,  the  approximation  appears  to 
conform  to  the  established  results  only  in  certain  cases. 
Some  possible  explanations  for  this  observed  deviation  will 
now  be  considered. 

Recall  that  the  integral  term  in  equation  (27)  was 
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dropped  (assumed  to  be  approximately  zero)  in  making  the 
Liouville-Green  approximation.  The  validity  of  this 
assumption  appears  to  be  linked  strongly  to  the  form  of 
P(z),  which  for  the  present  case  is 

m  n 

F(t)  =  [R  E  (m-n)  (3m  +  n+4)  ]  / 
a  b 

m+n  m+n+2 

[16a^b^(vel)  (t+ (R^-R^) /vel)  ].  (33) 

Hhen  the  exponents  m  and  n  are  equal,  F(t)=0,  and  by 
considering  equation  (27) ,  one  can  see  that  the 
approximation  is,  in  fact,  the  exact  solution.  This  is 
verified  trivially  by  Figures  3  and  4,  for  which  the 
exponents  are  equal.  It  should  be  noted  that  in  this 
situation,  the  model  reduces  to  the  case  of  a  constant  ratio 
of  attrition-rate  coefficients. 


When  the  exponent  m  is  greater  than  n,  F  (t)  will  always 
be  greater  than  zero,  and  from  (27)  it  follows  that  the 
exact  solution  is  greater  than  the  approximation.  This 
situation  is  illustrated  by  Figures  6  and  10.  The  converse 
also  holds.  Hhen  m  is  less  than  n,  F (t)  is  always  less  than 
zero,  and  the  approximation  is  greater  than  the  exact 
solution.  This  is  shown  in  Figures  5  and  9. 


Again  examining  (33) 


the  parameter  a  b  ,  which  is  the 

'00' 


product  of  the  opposing  forces'  fire  effectivenesses  at  r=0, 
appears  to  be  significant.  As  (a  b  |  increases,  F(t) 


decreases,  making  the  integral  term  in  (27)  smaller,  so  one 


would  expect  the  approximation  to  be  "closer"  to  the  exact 
solution.  This  did  occur  in  the  cases  considered  (see  for 
example  Figure  7) . 


B.  OFFSET  LINEAR  COEFFICIENTS 
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For  the  case  of  the  offset  linear  attrition-rate 

coefficients  (8)  ,  the  maximum  effective  range  of  the  Y  force 

weapon  system  (R  )  was  set  equal  to  1500  meters,  with  the 

a 

rest  of  the  parameter  values  equal  to  those  in  the  power 

coefficient  case.  The  maximum  effective  range  of  the  X  force 

weapon  system  (H  )  was  varied  throughout  the  range 

b 

considered  in  Ref  13.  The  results  are  shown  in  Figures  11 

through  15  where,  as  before,  each  figure  compares  the 
corresponding  solutions..  As  was  the  case  with  the  power 
coefficients,  the  solutions  are  seen  to  conform  only  in 
certain  cases. 

For  the  linear  coefficients,  it  was  shown  that  F  (z) 
could  be  written  in  terms  of  the  time  variable  t  as 

F(t)  =  E  R  (R  -R  )[12(t+(E  -R  ) /vel)  +  7  ( (R  -R  )/vel)  ]  / 
aBBa  aO  Ba 

3  3  3  ,  , 

£16a  b  vel  (t+ (R  -R  ) /vel)  (t+ (R  -fi  )/vel)  }.  (35) 

0  0  a  0  b  0 

When  the  maximum  effective  ranges  of  the  opposing  weapon 

systems  (R  and  R  )  are  equal,  F (t)  =0,  and  the  integral 
a  b 

term  in  (27)  is  zero.  This  again  reduces  to  the  constant 

attrition-rate-ratio  case  as  can  be  seen  from  (8)  with  A=0. 
Thus  one  would  again  expect  agreement  of  the  exact  and 
approximate  solutions.  This  was  found  to  be  true  as  shown 
by  Figure  11,  where  both  maximum  effective  ranges  equal  1500 
meters. 

As  was  the  case  for  the  power  coefficients,  the 
parameter  |a^b^|  appears  in  the  denominator  of  F(t),  hence 

the  two  solutions  should  be  closer  as  la  b  increases.  This 

'  0  O' 
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was  found  to  be  the  case,  as  illustrated  by  Figure  15. 


C.  ANALYTICAL  CONSIDERATIONS 

The  examinations  of  both  the  power  and  offset  linear 

attrition’-rate  coefficient  cases  have  centered  on  the 

parameter  F  (z)  in  the  solution  (27).  While  F  (z)  does  appear 

to  be  the  driving  factor  in  the  integral  term  in  (27)  (thus 

in  the  validity  of  the  approximation) ,  it  should  be  made 

clear  that  the  other  factors  in  the  integral  term  can  not  be 

ignored.  While  a  mathematical  analysis (for  example  of  the 

hyperbolic  sine  term)  would  be  formidible,  examples  can  be 

found  which  illustrate  that  F(z)  can  not  always  be 

considered  alone.  One  such  example  is  shown  in  Figure  8.  The 

fire  effectiveness  of  the  two  forces  at  r=0  were 

interchanged  (i.e.,  a  =.6  and  b  =.06  instead  of  vice  versa) 

0  0 

in  the  power  attrition-rate  coefficient  case.  From  the 

appropriate  form  of  F  (z) ,  equation  (33),  one  can  see  that 
this  change  has  no  effect  on  the  numerical  value  of  F(t) 
throughout  the  time  range.  Yet  the  significant  difference  in 
results  can  be  seen  in  Figure  8  when  compared  to  Figure  6. 
Currently,  there  is  no  convenient  methodology  to  analyze  the 
implications  of  such  discrepancies. 

Toward  this  end,  theoretical  analysis  in  the  area  of 
Giver's  "Error  Bounds  for  the  Liouville-Green  Approximation" 
[Ref  11 3  could  be  very  beneficial.  While  they  sound 
promising.  Giver's  results  can  not  be  directly  applied  to 
variable-coefficient  Lanchester-type  differential  equations. 
If  available  in  a  convenient  form,  such  error  bounds  could 
be  evaluated  for  specific  situations  to  investigate 
tolerable  error  limits  prior  to  invoking  the  approximation. 

From  the  above  analytic  considerations,  one  can  see 
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that  while  one  may  be  able  to  explain  some  of  the  error 
incurred,  this  is  small  comfort  to  the  analyst  attempting  to 
make  practical  use  of  the  Liouville-Green  approximation.  It 
has  been  shown  that  the  approximation  appears  to  yield 
credible  results  only  in  certain  limited  circumstances. 
Furthermore,  methodology  is  not  presently  available  to 
predict  the  error  which  will  be  incurred  in  a  specific 
situation,  except  for  the  trivial  case  in  which  the  integral 
term  in  (27)  goes  to  zero. 
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7III. 


SOHMAEY 


This  thesis  has  examined  the  ''adequacy''  of  the 
so-called  Liouville-Green  approximation  to  the  solution  of 
variable-coefficient  Lanchester-type  equations  of  modern 
warfare.  Although  this  approximation  only  involves 
"elementary"  functions  and  has  an  intuitively  appealing 
form,  it  was  unfortunately  found  that  this 
Liouville-Green-Lanchester  approximation  is  not  consistently 
reliable  for  , estimating  force  levels.  Specifically,  the 
approximation  was  applied  to  two  forms  of  attrition- rate 
coefficients:  (I)  power  attrition-rate  coefficients 
(modelling,  for  example,  the  same  maximum  effective  ranges 
for  the  opposing  weapon  systems) ,  and  (II)  offset  linear 
attrition-rate  coefficients  (modelling  different  maximum 
effective  ranges) .  It  was  seen  that  the  approximation  was 
made  by  dropping  the  integral  term  in  equation  (27) ,  hence 
exact  results  are  achieved  only  when  this  term  equals  zero. 
This  situation  was  found  to  occur  only  when  the  ratio  of 
attrition-rate  coefficients  was  constant.  For  the  power 
attrition-^rate  coefficients,  this  implies  that  the  exponents 
in  (7)  are  equal,  while  for  the  offset  linear  coefficients 
this  implies  the  offset  parameter  A  in  (8)  is  zero  (which 
occurs  when  the  maximum  effective  ranges  are  equal) . 
Analytic  considerations  showed  that  while  certain  initial 
parameter  combinations  produced  credible  results  from  the 
approximation,  there  was  no  methodology  available  to  predict 
the  error  incurred  in  the  general  case.  It  was  found  that, 
as  a  general  rule,  the  further  one  moves  from  the  constant 
coefficient  ratio  case,  the  larger  the  error  incurred  by 
invoicing  the  Liouville-Green  approximation. 
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TABLE  I. 


Properties  of  the  General  Lanchester  Functions  x^,  x^,  y^,  y^. 
1.  dXj^/dt  =  /k^/k^  a(t)y2 


2.  dx^/dt  =  \/\  a(t)y^ 


3.  dy^^/dt  =  /k^/k^  b(t)x2 

4.  dy^/dt  =  /k^/k^  b(t)Xj^ 

•5."^  x^(t)y^(t)  -  X2(t)y2(t)  =1  Vt 
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Table  II.  Properties  of  the  Power  Lanchester  Functions 


u  ,  V  ,  U  ,  and  V 
m,n  m,n  m,n  ir 


1.  du  /dt  =  /k^/k  [k  t”]V  (t) 

m,rL  b  a  a  iii,n 


2. 


dv  /dt  =  /k, /k  [k  t”^]U  (t) 

m,n  D  a  a  m,n 


3-  dU  /dt  =  /k  /k,  [k,  t'^lv^  ^(t) 

m^n  a  b  b  Di,n 

‘iV  /dt  =  /k  /k,  •  [Ic  t“]u  (t) 
111)11  a  b  b  ni)ii 


5. 


u  _(t)U  (t)  -  V  _(t)V  (t) 

111)11  ni)ii  111)11  01)11 


6. 


^(t=0)  =  U  ^(t=0)  =  1 

111)11  ni)ii 


7. 


V  (t=0)  =  V  (t=0)  =  0 
m,n  m,n 


8*  u  (t)  =  U  (t)  =  cosh  (/k  k, 

111)111  111)111  a  b 


9.  V  (t)  =  V  (t)  =  sinh  (/k  Ic 

TH)!!!  I11)I]1  a  b 


Properties  6  and  7  only  hold  for  m,n  >  -1. 


=  1  Vt 

^m+1^ (m+l) ) 

t"^V(nH-l)) 
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of  battle  is  denoted  as  R_  ®  1250  meters  and  (as  shown) 


X  FORCE 

tJ.OO  12- CC  16.00  2Q.Q; 


C3 


I 


"1 I 1 I I 

2.50  5.00  7.50  10.  CO  12.50 

fWMGE  two  )  ‘ 


Figure  3.  Power  coefficients  with  m=n=l. 

For  Figures  3  through  10,  R  =Rh  =  2000m .  , 
Xq=10,  yQ=30,  aQ=.06,  and  bf=.b  unless 
otherwise  specified.  '  ^ 
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X  FORCE 

-J.OC  li.OO  la.CO  ifi.OO  20.0! 


Figure  4.  Power  coefficients  with  m=n=2. 
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Figure  5.  Power  coefficients  with  m=2 ,  n=3. 
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Figure  6.  Power  coefficients  with  m=3,  n=2. 
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CO-yT  GC'?T  00 'F  OC't’ 

33ycj  v; 
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Figure  7.  Power  coefficients  with  m=3,  n=2,  a  =B  =.6. 

0  0 
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X  FORCE 

3.03  1C.  00  16  00  20-00 


Figure  8.  Power  co«;ff  i  ci  ents  with  m=3,  n  =  2,  a  =  .6,  B  =  .06 

0  0 
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pDRce; 


Figure  9.  Power  coefficients  with  m=2,  n=2.1. 
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x  foace 

d.oo  12.  oc  ir>.oo  za.oc 


Figure  10.  Power  coefficients  with  m=2.1,  n=2. 
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?t  FOflCE 

<1.00  8.00  12.00  16-00  20.00 


Figure  11.  Linear  coefficients  with  R  =  1500m. 

For  Figures  11  through  16,®  R^=1500m., 
Xo=10,  yQ  =  30,  aQ=.06,  and  t)Q=.6  unless 
otherwise  specified. 
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Figure  12.  Linear  coefficients  with  R  = 

B 


1 600m . 
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Figure  13.  Linear  coefficients  with  R  =  2000m. 
-  B 
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Figure  14.  Linear  coefficients  with  R  =  2500m. 
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Figure  15.  Linear  coefficients  with  R  =  3000m, 
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Figure  16.  Linear  coefficients  with 


Rg  =  2500m.  , 


a.  =  b  = 
0  0 
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